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1 Introduction, problem, results 

Anisotropic phenomena are plentiful in the natural sciences and technology. The list 
is long, we name just a few: crystal growth, phase transitions, image analysis, models 
of segmentation. We will concentrate on mathematical models involving operators like 
u t — divL(Vw), where L is a monotone map. Despite many year of research this area is 
full of challenging problems, the list of contributions is by no means exhaustive, [1], [3], 
[5], [9], [10], [12], [13], [14], [19], [20], [22], [31], [39]. 

In order to set a feasible goal we will concentrate on one dimensional version of these 
problems, thus we will study: 

u t -£L(u x ) = in(0,l)x(0,T)=/x(0,T)=:I T , 

u(0,t) = A, u(l,t) = B, 1 ' 1 

with a suitable initial datum u\ t=0 = u . In general, the boundary data are time depen- 
dent. Our basic assumption is 

L : R — > M is merely monotone increasing. 
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This is a minimal assumption to render (1.1) formally a parabolic problem. It also en- 
compasses a rather broad range of phenomena. It is enough to say that, we would like to 
address simultaneously the following equations, 

u t -u xx = 0, (1.2) 

u t - ^;(KI + u x ) = o, (1.3) 

u t - -7- (sgn (u x )) = 0. (1.4) 
ax 

We see that mere monotonicity of L means that we admit jumps and flat regions of L, 
which represent singular and degenerate parts of the operator. 

We notice that in equation (1.3) the diffusion turns on only for u x > 0, i.e. for areas 
where u is increasing. There is no evolution for decreasing parts. This seems to be the 
simplest model with anisotropic character. On the other hand, in equation (1.4) the 
diffusion is present only in the zone, where u x = 0. We know (cf. [17]) that L should be 
understood as a maximal monotone multivalued operator, e.g. L(p) = sgn (p) is given by 



sgn (p) 



We want to establish in this paper a general existence result and at the same time the 
maximal possible regularity permitted by the generality of our assumptions on L. 

In principle, one can claim that existence and uniqueness is well-known. Our maximal 
monotone operator L is a subdifferential of a convex function W, L = dW. Next, we set 

T , . _ J fj W(u x ) dx if the integral is finite, , , 

+00 elsewhere. 

Thus, equation (1.1) is a gradient flow of J in L p and we can apply the abstract nonlinear 
semigroup theory, see [7], [17]. This approach has been implemented by many authors, 
including [1], [3], [21]. We can call the solution, we obtain in this way, variational, because 
constructing the Yosida approximation requires solving a variational problem. 

However, discussing regularity with this tool is cumbersome, because it requires de- 
termining the domain of the operator 

A:D(A) C L 2 (I)^L 2 (I), 

given by 

d 

A(u) = —L(u x ). 

Of course, A has to be properly understood. 

An additional complication comes from the fact that the boundary data may be time 
dependent. In any case, there is no common description of D(A). It suffices to consider 
in (1.1) A = = B, and to mention that if L(p) = p, (i.e. we have (1.2) in mind), 
then D{A) = H 2 (I) R Hq. On the other hand, when L(p) = sgn(p), then D(A)<^{u G 
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AC : u\qi = 0, u x £ BV}. In addition, we know that D{A) does not contain all 
smooth functions, e.g. we know that x 2 ^ D(A) (cf. [37]). As a matter of fact a 
complete description of D(A) for that L is not available yet. However, it is available for 
the multidimensional version of (1.4), see e.g. [4]. It is based on Anzellotti's formula 
for integration by parts, [6]. However, a direct characterization of this set for the one- 
dimensional problem seems to be missing, a partial result in this direction, for Neumann 
boundary data, is in [30]. 

Another aspect of the same problem is shown in the papers by Meyer (see [32]) and 
Mucha (see [34]). These authors showed that for A the total variation operator, i.e. 
Au = div and / ^ but small, the only solution to Au = f with the homogeneous 
Dirichlet boundary condition is u = 0. 

Therefore, it is not surprising that the complete discussion of regular of solutions is 
beyond the scope of this paper. 

It is worth noticing a development of the viscosity theory for equations like 

u t = a(u x ) (W p (u x ) x + a) , 

where W is merely convex. This line of research was initiated by Giga and Giga, see [23]. 
The point is, this theory is capable of handling faceted solutions too, but facets have to 
be of positive length. In [24] a comparison principle was proved for such solutions and a 
number of examples was presented. However, there is no general existence theory, yet. 
Before we discuss the regularity of solutions, we state the basic existence result. 

Theorem 1 Let us assume that L(-) is monotone increasing, A,B £ W^OjT), u £ 
Li(0, 1), such that uo tX £ BV(0, 1), then there exists a unique solution u to problem (1.1), 
i.e. 

u x £ ^(0, T; BV[0, 1]), u(0, t) = A and u(l, t) = B 
and there exists Q £ L 2 ((0, 1) x (0,T) such that following integral identity is valid 

K, 0)l 2 (o,i) + ( n » 0x)l 2 (o,i) =0 for a.e. t>0 (1.6) 

and for all <p £ C °°((0, 1) x (0, T)) and 

fl(x, t) £ L o u x (x, t), (1.7) 

where the last term is treated as a composition of multivalued functions. 

In fact, we shall see that the space derivative d x u(-,t) £ i?V[0, 1] for t > is the 
highest possible regularity, permitted in such generality. We are able to show that only 
because we consider a one dimensional problem. Of course, such a general result need not 
be optimal for all choices of L. 

In order to describe the regularity of solutions we will use here the notion introduced 
in [37], [36], and developed in [26]. Namely, we introduced the notion of almost classical 
solutions. Roughly speaking, u, a suitably defined weak solution is almost classical if 
it satisfies the equation pointwise, except a very small set, see Definition 3. In [37] we 
showed that the equation, we studied there, is satisfied except finitely many points. We 
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will illustrate this notion assuming the initial conditions are convex. This restriction is 
for the sake of simplicity only, originally in [37], [26] we do not impose such restrictions. 

The difficulty, related with this notion, is that it is not quite pointwise. In reality we 
have to give a meaning to the composition of L and u x , when both of them are discon- 
tinuous. It turns out that it is convenient to interpret these two functions as multivalued 
operators. While is it quite easy to interpret L as a maximal monotone operator, more 
work is necessary if u x G BV . In this case u is Lipschitz continuous and we may use 
the notion of Clarke differential of u to define set- valued u x . The details on the Clarke 
differential can be found in the Appendix. 

Of course, composition of multivalued monotone operators in general need not be 
monotone. However, we solve this issue by introducing Lou x , which is a special selection 
of Lou x , see Definition 2. The main advantage is that indeed Lou x is a monotone operator, 
and on facets it is even single valued. The drawback of this definition is that it is still 
essentially one-dimensional. 

The existence result, Theorem 1, involves passing to the limit in L e (u e x ), where quan- 
tities indexed by e are related to an approximate system, when both limiting functions L 
and u x may possibly be discontinuous. Thus, it is worthwhile to look more closely at the 
limit of L e (u e x ). The question is to which extent we can call Q, the limit of L e (u x ), the 
composition L o u x or rather flclo u x . In principle there is a lot of choice for Q, but by 
uniqueness theorem the freedom is rather limited. Finally, it turns out that we can give 
a formula for Q that based on an explicitly defined composition Lou x . 

A few aspects of (1.1) have been well studied. The most prominent example is the 
case of uniformly parabolic equation, when L is smooth and j^(p) > e > 0, see [29]. 
The degenerate problems have been studied in [1] and from a different perspective in [8]. 
The TV-flow is another extreme case. It has been well studied due to its connection to 
image processing. There, frequently, the domain is rectangular and the solution satisfies 
the homogeneous Neumann boundary conditions, see [3], or it is studied in the whole W 1 , 
[2], [11]. However, the case of Dirichlet data has been also studied, see [4]. The studies 
motivated by crystal growth are by [21], [25], [27], the anisotropic version of the total 
variation flow (1.4) was performed in [33]. The TV flow with constraints was analyzed in 
[22]. Finally, the one dimensional version of the total variation flow was studied in [26] 
and [16]. Some basic questions have been studied for a system with L(jp) = p + sgnp 
in [35]. However that paper showed just the direction of our research program, not a 
complete general theory. The present paper provides the missing existence result. 

The complete analysis of (1.1), taking into consideration all possible types of behavior 
of L is beyond the scope of the paper. However, we may relatively easily discuss the case 
of Mo being a convex function. The result below guarantees that convexity is preserved. 

Theorem 2 Let the assumptions of Theorem 1 are fulfilled. In addition. A, B are time 
independent and Uq is a convex function. Then for each t G (0,T) the solution u(-,t) is a 
convex function, too. 

By Theorem 2, we may restrict our attention to convex data uq. Since at each time 
instance u x (-,t) is a monotone function, we can extended to a maximal monotone graph. 
We define such composition Lou x of maximal monotone graphs L and u x so that Lou x is 
a monotone operator. 



4 



Now we want to distinguish two sets for / and L via a definition below. 
Definition 1 Let u be a convex function over I, then we set 



D f (u T ) 



{\J[**M ■ Ux\[a k ,b k ] = °k for a constant 9 k and a k < 6 fc j. (1.8) 



'f\u x 

k 

Let L be an increasing function over R, then we set 

S L = {u : L(uj) = [a, b] with a < b}. (1.9) 

The set Df(u x ) describes the flat parts of u x and Sl jump points of L. 

Now, we are prepared to introduce the definition of our special composition o. 

Definition 2 Let L be a multivalued operator given by (1.6) andu x G BV(I) be a deriva- 
tive of a convex function. Then, the multivalued operator 

Lo Ux (1.10) 

is defined as follows: 

1. Ifu x (x ) i S L , then 

Lou x (x ) = L o u x (x ). (1.11) 

In particular, ifu x {xo) = [p,q], then Lou x {xq) = [L(p),L(q)]. 

2. If u x (xq) = 9 G Sl and L{6) = [a,b], then: 

2.1. ifu-\6) = [£_,£+] with < £+ and £_ ^ 0, £+ ^ 1, then 

Lo Ux {s) = x k {s- £+) + y k (s - £_) for s G (£_, f+) (1.12) 

with 

a b 
Xk = 1 T and Vk = 1 T'i K 1 - 1 ^) 

2.2. ifu~ l (0) = {x Q } ; then 

Ldu x {x ) = [a,b); (1.14) 

2.3. ifu x 1 {6) = [0,f+] £+ G [0, 1) (or [£_, 1] wift G (0, tfien we put 

Lou x {s) = b for se [0,£ + ] (Lou x {s) = a for s G 1]); (1.15) 

S.^. if u~ 1 (9) = [0, 1], £/ien we put 

Lo Ux (s) = (3 for se[0,l]. (1.16) 

By definition, Lo-u^. =: <j is an element of a composition of two multifunctions, L o u x . 
In the examples below, we frequently specify a selection a, which is absolutely continuous. 
Thus by the uniqueness theorem, it is sufficient that the equation (1.1) is satisfied point- 
wise a.e. Observe that the composition Lou x is maximal monotone when u is a convex 
function. 

Finally, we define almost classical solutions. 
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Definition 3 Let us suppose that u is a solution given by Theorem 1. We say that u is 
an almost classical solution to system (1.1) iff the following identity holds 

u t -±Lou x = in Ix(0,T))\E, (1.17) 

where 'H 1 {E) is finite. 

The examples we constructed so far (see [37]) show that we cannot expect that E be 
smaller than stated above. 

The definition of the almost classical solutions, introduced in [26], is based on the 
composition 5. We refer to this paper for details. Here, we present the definition which 
covers the needs of our analysis, namely the case of u being a solution with convex initial 
data. We shall emphasize that the set of irregular points is at most countable for t > 0, 
but if L suffers just a finite number of jumps, or more precisely, jumps are isolated, then 
the number of irregular points is finite for each t > 0. 

The final result clarifies the importance of solutions to (1.1) with convex data. 

Theorem 3 Solutions with convex initial data coming from Theorem 1 are indeed almost 
classical. 

The rest of the paper is devoted to the proofs of these results. In the next section 
Theorem 1 is shown, in Section 3 we prove results for convex solutions: Theorems 2 and 3. 
Some parts of proofs are stated very precisely, although they belong to the classical/well 
known theory. But they are important to explain the need of redefining the meaning of 
solutions as well as to initiate the revision of regularity of solutions. 

2 Existence and uniqueness of variational solutions 

In this section we prove Theorem 1. We proceed in a standard way. The monotone 
operator L, will be substituted by a strictly monotone, single valued smooth function 

L*(p) = (L*n%p) + ep, (2.1) 

where ir e is a smooth non-negative approximation of the Dirac delta. 

Now, we construct an approximating systems. For a given e > we examine 

uf = m (0,l)x(0,T), 

u e (0,t) = A e , u(l,t)=B e for te(0,T), (2.2) 

u e \ t=0 = u e Q on (0,1). 

Data A e ,B e ,UQ are so regularized that they are smooth functions on their domains and 
the consistency conditions are satisfied. 

Due to the classical theory, see [29], system (2.2) admits a unique smooth solution, 
at least locally in time. Hence, we concentrate our attention on finding suitable a priori 
bounds, which allow us to pass to the limit with e — > 0. Existence and uniqueness will be 
established on a given interval [0, T]. 
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Let us formulate the basic energy estimates guarantying the existence solution to the 
approximating system, stipulated by Theorem 1. For this purpose we reduce the given 
problem to one with homogeneous boundary conditions. We set 

u e - [(B e (t) - A e (t))x + A e (t)} =: v e . (2.3) 

Let us subtract [(B e (t) - A e {t))x + A e {t)] t - £L e ([{B e (t) - A e (t))x + A e {t)] x ) from both 
sides of (2.2). Then, we have 

vf - £[L*(u%) - L<(&(t) - A'(t))] = 

-[((B*y(t) - (A*)'(t))z + (A«)'(t)] + £L<{B<(t) - A<(t)) in (0, 1) x (0, T), 

v e (0,t) = 0, u e (l,t)=0 for te(0,T), ^'^ 

v% =0 = ul - [(B £ (0) - A e (0))x + A%0)] on (0, 1). 

Next, we test the equation (2.4) with v e . 

~ Jyfdx + j[L*{u x ) - U{B\t) - A\t))\[u% - mt) ~ A\t))]dx = 

- [(m(t) - {A*)\t))x + (^)'(i)] + ^L\B\t) - A<(t)j\ v*dx. (2.5) 

We notice that for fixed, positive e the value of ^L e (S e (t) — A € (t)) is clearly zero, because 
L e is single valued. 

Identity (2.5), combined with monotonicity of L e , leads to the following estimate 



d 1 
dt2 



iv e ydx < / [-[HB e )'it) - iA e )'it))x + iA e ) f it)]v e dx. 



Taking into account Gronwall inequality yields a desired e-independent estimate, 
IKWIlLm < (lKlk(/) + / 

0<t<T 



sup Mt)\\Un ^ (lKlk ( i)+ [ T \\((A*)'(t)-(B*y(t))x- (AJ(t)\\l 2(I)) dt 

D<t<T V JO 



(2.6) 



Let us notice that if A and B are time independent, then we will obtain a better 
estimate in place of (2.6), i.e. 

SU P IkllLm < KIU 2 (/)- ( 2 - 7 ) 

0<t<T 

Now, we would like to find a bound on the second space derivative of u e . After having 
differentiated (2.2) twice with respect to x, we will have, 

C-^(^ £ KK,) = in It, 

£ p L*(ul)u xx (0,t) = A$, £ p L*(u x )u xx (l,t) = B\ for t e (0,T), (2.8) 

Kx\t=o = u € on (0,1). 
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The boundary data are obtained directly from the system (2.2) by the analysis of the 
boundary values of u\. By the analysis of integrals over set {u xx > 0} and {u e xx < 0} we 
shall momentarily show a form of the maximum principle 

ll M LIUoo(0,T;Li(0,l)) < IhcUJUlCO.l)- ( 2 - 9 ) 

Let us fix e. Despite smoothness of u e we cannot claim that the set {u xx = 0} is 
regular. However, Sard theorem gives us a sequence a k — > such that, indeed, each 
{u xx = a k } is regular. 

Let us fix a k . In order to examine set {u e xx (x) > a k } we look at (2.8) written in the 
following form, 

Kx ~ *k)t - i(^ f (^)(«» - (Tk)) = eh&iL<(ul) in (0, 1) x (0,T), 
|L e (u«H*(0,*)=A e , T p L e (K)u xx (ht) = for t e (0,T), (2.10) 

u xx \t=o = u e 0:XX on (0,1). 

Let us integrate equation (2.10i) over {u e (x) xx > a k }. Noticing that 

d f 
dtJi 



(u e (x) xx - a k ) + dx = I (u e (x) xx - a k ) t dx 



yields after integration by parts, that 

d f f d d 

-77 / {u xx -a k )dx- I n d{u c xx _ ak} (x) — (—L e (u x )u xx -a k )dfi 

at J{u% x >a k } Jd{u xx >a k } dX dp 

=°*L^M L ' i<)ix - (2 - n) 

Here, we denote the atomic measure on d{u xx > a k } by dp,, the term rid{ u e >ak \(x) is 
the normal vector, here it is ±1, indicating the orientation. The boundary integral is 
well-defined due to regularity of d{u xx > a k }. Moreover, we can bound it too. Indeed, if 
for xo G d{u xx > a k } and ng{ u e > ak }(xo) = — 1, then for x G (xo,xo + 5) we are ensured 
i^(ul)(ul x -a k )>0, hence 



that 4zL e {u%){u xx -a k ) > 0, hence £{±L e {u e x )(u xx - a k )){x ) > 0. By the same token, if 



n d {u xx >a k }{xo) 



then for x G (x — 5, x ) we are ensured that dL e (u x )(u xx — a k ) > 0. As a result, 



:(dL e (u x ) (u xx — cr k ))(x ) < 0. Thus, we conclude that 



f d d 

/ n a{u%x>ak] (x) — (—L e (u x )(u xx -a k ))dji>®. (2.12) 

Jd\u%„.>o h \ ax a P 



'9{u% x >a k } 

This is why we arrive at 

-T r A2 
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We recall that sequence Ok was chosen for fixed e > 0. We are allowed to pass to the limit 
with cjfc. Since the solutions as smooth, the integral in the right-hand-side (rhs) is well 
defined and finite, so after taking the limit — > we get 



sup / (u xx ) + dx < / (u 0:XX ) + dx. (2.14) 

0<t<TJl J I 

Repeating this argument for {u xx )_ = max{— u^O} yields 

sup / \u xx \dx < / \u e \dx. (2.15) 

0<t<TJl J I 

We aim at using the lower semiconituity of the BV semi-norm with respect to L\ 
convergence. For this purpose we have to establish that u £ x tends to u x as e — > 0. In order 
to capture this we need a bound on u\. We go back to the equation and together with 
the information implied by (2.2) we obtain 

KIIl^T;^/)) < l|£ e (OIU»(0,T;Loo(J)) < M. (2.16) 

The rhs of (2.16) is uniformly estimated, because by (2.15) u x is pointwisely bounded and 
from the definition of L e we deduce that L e ([— 1 — Hw^Hloo, 1 + IKllUooD is a bounded set 
too. 

On the way, we would like to obtain some information in the Holder class, which is 
of independent interest. For this purpose, we have to modify a bit the sequence. Take 
the extension operator onto the whole line M, such that the function will be compactly 
supported and the considered regularity will be preserved too. Namely, we take U e = Eu £ . 
Next, we define 

V e = (l-dl)- 1/2 U e . (2.17) 
Estimates (2.15) and (2.16) combined yield the following bound for any finite p > 1, 

V; e L p (0,T;L p (R)), V£ G L P (0,T; L p (R)), i.e. V £ G W^(R x (0,T)). (2.18) 

The embedding yields W^ 1 C C 2a ' a (R x (0, T)) with a = 1 - ^ > \, provided that for 
p > 3, [15, Chap. XVIII]. Going back to u e we get 

u e e G 2a-l,a-\ ^ ^ 

This clarifies the strong/pointwise convergence of u itself, but also we need information 
about the strong convergence of u%. From (2.15) we deduce an extra regularity due to 
u% G L p (0, T; Wp /p (I)). We observe that ^(0, T; W?(R)) C L p (0, T; W p 1+1/p (R)), because 
|(1 - i) + = i ; [15 ; Chap.XVIII]. Thus, due to the definition of V e we obtain that 

V £ G L p (0, T; H/ p 2+1/p ' 1 (M)) and V* xt G L p (0, T; W p - 2 (R)) for 1< p < oo. (2.20) 
We can restate (2.20) in terms of U to obtain 

V% G L P (0,T; ^(1)) and U xt G L p (0,T; W~^ P (R)). (2.21) 
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Since W^ p C L p C W_ 1 and the Aubin-Lion theorem applied to u e defined over / yields 
(up to a subsequence — > 0) 

< -» <. strongly in L p (0, T; L P (J)). (2.22) 

Passing to the limit e — > 0, while using (2.15) and the lower semiconituity of the BV 
seminorm yields, see [40], 

ess sup \\u x (-,t)\\ BV(I) < \\u 0jX \\bv(I)- (2.23) 

0<t<T 

Due to (2.19) we control the solution at the boundary and 

u € \ai -> u\di strongly in C[0,T], 

and of course the boundary conditions (1.1)2 ar e fulfilled. 

Once we established (2.23) we want to use it to study u\ and the composition L € (u e x ). 
Our goal is to show that for almost all t, L e (u e x (-,t)) has a pointwise limit denoted by 
Q(-,t). The problem is lack of continuity of L. 

Let us first set g e (x,t) = (B e (t) - A £ (t))x + A e (t), then v e satisfies 

vt - ^L\ul) = -g\. (2.24) 

Now, we multiply (2.24) by v\ and integrate over Q^- After integration by parts the 
second term on the l.h.s. and using that v e \gi = we obtain, 



T r rT 

e\2 



[(vty + L*(uiyjdxdt = - g\v\dxdt. (2.25) 

I Jo J I Jo 

Let us now introduce W € ClS db primitive function of Z/, 

W e (p) = f L e (r)dr + C, 



where constant C is so chosen that for allp G [— 1 — H^Hl^, 1+|| m xIUoo] we have W(p) > 0. 
Monotonicity of L e implies convexity of W e . Hence, due to the definition of v e we obtain 

2~ / / ^ 2dxdt + f we ( u x) dxdt 

< J^W e (u , x )dxdt + ^JJ o (9tT dxdt + J J L e (u e x )g € t dxdt. 
The definition of v e yields, 

I [ (u\) 2 dxdt< ! ! (glfdxdt + A ! W e (u ^ x )dxdt+ f [ (L e (u%)fdxdt. 
J 1 Jo J 1 Jo J I J I Jo 

Boundedness of ||i4lkoo implied by (2.23) yields that the last term on the r.h.s. above is 
finite. Thus, we have shown that 

KXlk(n T ) < C(\\A*\\ wHOjT) , \\B*\\ w i (0tT) ,L, IKsllw), (2-26) 
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where C is independent of e. We now conclude that 

u e — ^ u in W^(I T ). 

At the same time, by (2.16) L e (u e x ) is bounded in L 2 (It), hence after selecting a subse- 
quence 

L e (u e x )^n in L 2 (I T ). (2.27) 
We will show that the convergence is better. 

Lemma 2.1 The limit functions u and Q fulfill (1-6), i.e. for almost allt>0 the identity 
u t = fir holds in the sense of distributions on (0, 1). 

Proof. We notice that for a fixed ip G C^(I) and any t Q G (0, T), r > expression 

J I u\ipdxdt 

Jt -T JI 

is a continuous functional over W%{It)- Then, weak convergence of u\ implies that 

to+r r rto+r r 

/ u\ipdxdt — >■ / / u t (pdxdt. 
t -T Ji Jt -r ji 

Using the equation (2.2) and (2.27) we shall see that 

fto+r r rto+T r rto+T r rto+r r 

/ / = / / m Ux )) x <p = - / u[u x )^ x / / thp x . 

Jto-T J I Jto-T JI Jto-T JI Jto-T JI 

Thus, for all t G (0, T), r > and </? G C$°(I) we have shown that 

rto+r r rto+r r 

- tt(p x = / u t <p. (2.28) 

7t -T JI Jt ~T JI 

Since is fixed, then Lebesgue differentiation theorem yields 

— J Q(x,t)ip x (x) dx = J u t (x,t)ip(x) dx, for a.e. t G (0, T). (2.29) 

The choice of t depends on (p. We repeat the argument from [37, Theorem 2.1, page 2292] 
to claim that (2.29) holds independently of <p on a set of full measure Q C (0, T). 

Identity (2.29) means that t) has a weak derivative and fi x -(-, t) = u t (-, t) for almost 
all t. In addition, 

||^IU 2 (0,T;H 1 (/)) ^ |Klk 2 (ir)- 

Using the argument as in [37, Proposition 3.2] we come to the following conclusion. 

Lemma 2.2 The limiting function Q is a section of the subdifferential dJ at u (functional 
J is defined in (1.5)), i.e. 

J W(u x + h x ) - W(u x )dx > J Qh x dx. (2.30) 
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Proof. Because of the convexity of functions W e , we have the inequalities, 

J{W\u% + h x ) - W'(uD) dx>£ ^-W\u%))h x dx, 
where h G C£°(I). We integrate both sides with respect to t, thus we obtain 

Y T [ t+T J( we K + M - W '( U D) dxdt > y t f t +T I ip we «)) h * dxdt > ( 2 - 31 ^ 



We pass to the limit on the r.h.s. of (2.31). If we take into account (2.22), then we 
can apply Krasnosielskii theorem, see [28], to the l.h.s. of (2.22) while keeping in mind 
uniform boundedness of ||i4lUoo given by (2.23). This results in the following inequality, 

rt+r i> rt+r i> 

27 J J ^ W{Ux + kx) ~ W{Ux) ^ dxdt ~ 27 J J Qkx dxdt 

Using the argument as before we claim that (2.30) holds for a.e. t G (0, T). 
We show more, 

Lemma 2.3 If (2.30) holds for a fixed t and all h G H%(I), then Q(x, t) G dW o u x (x, t) 
for all x, except at most countably many. In other words Q(x,t) G ou x (x, t) a.e. 

Proof. Since u x (-,t) G BV(I), then u x is continuous except for at most countably many 
points. Let us suppose that u x (-,t ) is continuous at (x ,t ) and p G M is arbitrary. We 
take h € G Cg°(I) such that 

supp/i e C (xq — e,x + e), h x (x ) = p, \h%\< M, 

where M is independent of e. 

Now, we notice that for a function g continuous at xq we have ^ J x °_* g(s) ds — > g(xo). 
This observation and (2.30) combined yield 

W(u x (x , t ) +p)~ W(u x (x , t )) > fi(zo, to)p, a.e. x , t , 

as desired. □ 

Proposition 2.1 Variational solutions to (1.6) are unique. 

Proof. In order to show the uniqueness of our solutions it is enough to work in the 
L 2 -framework. We first notice that the difference of two solutions is an admissible test 
function. This is so, because the difference vanishes on the parabolic boundary of It- 
Take two solutions to (1.6), say u\ and u 2 , then we have 

(«i,t - u 2 , t , 0) + (fii - n 2 , <f>x) = in / x (0, T) (2.32) 

and u\(0, t) — 112(0, t) = and ui(l, t) — u 2 (l, t) = for t G (0, T). Thus, we insert it into 
(2.32), 

-f, 2 ( M i _ M 2, ui - u 2 ) + (fii - fi 2 , u ljX - u 2 , x ) = 0. (2.33) 

By Lemma 2.3, we deduce that (Tlx — Q 2 ,u\^ x — u 2 , x ) > 0, since the initial data are the 

— I 

dt 



same, so (2.33) yields 4||^i — u 2 \\ 2 < 0, i.e. u\ = u 2 on the full interval of the existence 



□ 
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3 Convexity of solutions 



In the previous section we showed existence and uniqueness of solutions to (1.1). Here, 
we prove results for convex solutions, i.e. Theorems 2 and 3. We proceed as in the proof 
of (2.9). 

Proof of Theorem 2. We aim at proving that for t > 0, solution u(-,t) is a convex 
function. Since u is obtained as a limit of approximating solutions u e to (2.2), it is 
sufficient to show that for alH > and e > functions u e (-,t) are convex. 

Let us consider system (2.2). Since u e at the boundary is fixed and we assumed that 
A, B are time independent, then u\ = at dl, so ^L e (u x )u xx = at dl, too. As a result, 
we find that u xx = at dl. Making use of the fact that u e are smooth, we are allowed to 
differentiate twice eq. (2.2). This yields 

<**-fe[^ e KKJ = m /x(0,T), 

u e xx = at <9/x(0,T), (3.1) 

u xx \ t =o = ul xx > at /. 

Although u e is smooth, we cannot claim that the sets {(x,t) e / x [0,T] : u e xx (x,t) = 
0} are regular. However, the Sard theorem guarantees us that there exists a sequence 
5k — > 0~, such that sets {(x,t) G I x [0,T] : u e xx (x,t) = 5k} are regular submanifolds, 
where 5k < 0. We follow the notation /_ = max{— /, 0}. Next, we restate (3.1) as follows 

- - - «1 = S *&4 L ' K) - (3 ' 2) 

Next, we integrate (3.2) over the set Ak = {(x, t) & I X [0,T] : u xx (x,t) < 5k} 

QA i< x ~ h)n t da - ^[^K)ul x - 5 k )]n x da = 5 k £ ^^K). 

Here, (n x ,n t ) is the outer normal vector to A\.. 

Since the set dA k is regular and the function u xx equals 5k on dA k \ I x {0,T}, then 
we § et JdA k ( u *z ~ 6 k)n t da = 0, moreover, j dAk £[-^L e {u e x )u xx - 5 k )]n x da < 0, because 
of ^L e {u x ) > and dl n (0A k Dl x {t}) = for t E (0,T). The initial datum u e 0xx is 
non- negative and 5k < 0. Thus, we obtain 

J{u xx -5 k )_dx<\5 k \ £ J\^-{± L \u x ))\dxdt. (3.3) 

After letting k —)■ oo ( 5 k — > we get 

y«J_dx = 0. Thus < x (x, t) > for t > 0. (3.4) 

Next, after passing with e to 0, we complete the proof of Theorem 2. 
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Proof of Theorem 3. We proceed in a few steps. We know that ut G Li{Jt) hence, for 
almost all t, with respect to one-dimensional Lebesgue measure, we have u(-,t) G L 2 (I). 
We have also established that for almost all t the equality 

u t = £l x (3.5) 

holds in the L 2 sense. In particular Q G H l (I). 
Here comes the main observation. 

Lemma 3.1 Let us suppose that Uq is convex and u is a corresponding weak solution to 
(1.1), then for almost all t function f2(-,£) satisfies Q(x,t) = Lou x (x,t). Hence, Q is 
differentiable except at most countably many points. 

Proof. We have to consider the cases specified in the definition of o. The crucial one is 
when for a number c we have L(c) = [a, b] with a < b. From Theorem 2, we know that 

{x G / : u x (x,t) = c,t > 0} = [£_(*),£+(<)] (3.6) 

where £_(£) < since u(-,t) is convex. 

The monotonicity of L(-) and u x (-,t) implies that we are able to find two sequences 

i X n (*)}) i X n (*)} SUCn tnat 

-> r(t)" and xt -> e + (t) + (3.7) 

as well as 

L(u x (x~(t),t)) ->■ a and L(u x (x+(t),t)) ->■ 6. (3.8) 

Namely, we require that points a^(t) are regular of function u x (-,t) and u x (x^(t), t) are 
regular of £(■), too. Then, we look at 

= t) on x+(t)]. (3.9) 

The integration over the marked set yields 

= L(u x {x+{t),t)) - L(u x {x-{t),t)). (3.10) 

x„ (t) 

The form of the rhs follows from the choice of suitable regular points. Since we are allowed 
to pass to the limit with n, we get 

C+(t) 

n^x = b — a. (3-11) 

£-(*) 

On the other hand, function u restricted to the set [£_(£),£+(£)] is linear, thus it has 
the following form 

u I [?-(*)>£+(*)] = cx + ( 3 - 12 ) 
Note that the facet is prescribed by the value c, so the slope is fixed in time, however the 
shift given by a(t) may depend on time. Of course, this holds only on the considered set. 
Hence, (3.12) follows 

"t I [€- (*M+(t)] = «*(*)■ ( 3 - 13 ) 
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Combining this with (3.11), we get the following relation 

at(£+(<) -£-(*)) = &-<*■ (3.14) 
By Definition 2, we observe that 

Hence, comparing (3.5), (3.14) and (3.15) we found 

ut = ^,Lou x at g_(t),e+(t)]. (3.16) 

On the remaining parts of the domain we have the regular case. By Lemma 2.3, we 
are ensured that fl G L o u x , and since by (3.8) L(u x (x~ (t) , t)) = Q(x~(t),t) — > a and 
L(u x (x+(t),t)) = £l(x+(t),t) —> b, and Lou x (£~(t),t) = a and Lou x (£ + (t),t) = b. So, we 
get 

£l(x, t) = Lou x as i? 1 (0, 1) functions. (3-17) 

We should stress that simplicity of the consideration follows from the convexity of function 
u(-,t). The proof is done. □ 



Appendix 

In (1.7) we compose two multivalued operators L with u x . Here we explain how we 
can treat BV functions as multi-valued operators. This is easy for functions which are 
derivatives, u x G BV{I). It is very useful in the regularity study of solution to (1.1). 
If u and u x belong to BV(I), then u is Lipschitz continuous. Hence, and exist 
everywhere and they differ on at most countable set. Thus, we may set 

d x u(s) = {ml + (1 - t>+ : t G [0, 1]}. (3.18) 

Under our assumptions on u, the set d x u(x) is the Clarke differential of u and equality 
holds in (3.18) due to [18, Section 2, Ex. 1]. If u is convex, then d x u is the well- 
known subdifferential of u. As a result, if u x G BV, then for each xq G (0,1), we 
have d x u(x ) = \\xm x _^ x - u x (x), lim a ._ +a .+ w x (a;)] or , where [a, 6] or = [a, b] for a < b and 
[a, b] or = [b, a] for b > a. 
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